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An aeroelastic analysis of the behavior of an entry vehicle utilizing an attached inflatable 
aerodynamic decelerator during supersonic flight is presented. The analysis consists of a 
planar, four degree of freedom simulation. The aeroshell and the IAD are assumed to be 
separate, rigid bodies connected with a spring-damper at an interface point constraining the 
relative motion of the two bodies. Aerodynamic forces and moments are modeled using 
modified Newtonian aerodynamics. The analysis includes the contribution of static 
aerodynamic forces and moments as well as pitch damping. Two cases are considered in the 
analysis: constant velocity flight and planar free flight. For the constant velocity and free 
flight cases with neutral pitch damping, configurations with highly-stiff interfaces exhibit 
statically stable but dynamically unstable aeroshell angle of attack. Moderately stiff 
interfaces exhibit static and dynamic stability of aeroshell angle of attack due to damping 
induced by the pitch angle rate lag between the aeroshell and IAD. For the free-flight case, 
low values of both the interface stiffness and damping cause divergence of the aeroshell 
angle of attack due to the offset of the IAD drag force with respect to the aeroshell center of 
mass. The presence of dynamic aerodynamic moments was found to influence the stability 
characteristics of the vehicle. The effect of gravity on the aeroshell angle of attack stability 
characteristics was determined to be negligible for the cases investigated. 
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Nomenclature 

= initial amplitude 
= decay constant 

= aerodynamic axial force coefficient 

= aerodynamic pitching moment coefficient slope about attachment point, dC m fda 
= aerodynamic pitch damping coefficient about attachment point, dC m / d(f)d /2V) 
= aerodynamic normal force coefficient slope, dC N jda 

= vehicle diameter 
= aerodynamic axial force 
= aerodynamic normal force 
= acceleration due to gravity 
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= angle of attack 
= interface damping 
= non-dimensional interface damping 
= interface stiffness 
= non-dimensional interface stiffness 
= angular frequency 
= phase angle 
= atmospheric density 
= surface atmospheric density 
= pitch angle 

= pitch difference between aeroshell and 1AD 


= superscript indication of IAD quantity 


CL = Centerline 

IAD = Inflatable Aerodynamic Decelerator 

LH = Local Horizon 


I. Introduction 

The recent resurgence of research activities related to the development of inflatable aerodynamic decelerators 
(IAD) aims to create an enabling technology for global Martian surface access and the landing of large payloads, 
especially for vehicles with high ballistic coefficients. Although the beneficial drag characteristics of IADs are well- 
established, 1 the question of stability for various shapes has yet to be answered satisfactorily. Knowledge of the 
stability of an IAD is important for ensuring the vehicle remains within a tolerable range of rotation rates during 
flight. This study is a first-order analysis of the stability characteristics of a flexibly-connected aeroshell and IAD. 

As the interface between the aeroshell and the IAD is flexible, an aeroelastic multibody simulation is required to 
capture the effect this flexibility has on the free flight stability characteristics of the system. This analysis is an 
aeroelastic study of a coupled, flexible aeroshell-IAD system in supersonic free flight. The development of the four 
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degree of freedom equations of motion is described along with the aerodynamics involved. Non-dimensional 
parameters are introduced to assist in the evaluation of the results. The simulation is validated in the rigid body limit 
with previously published rigid body simulations and with angle of attack amplitude growth theory. 2 Two flight 
regimes are studied using the simulation. The first is a constant velocity simulation where the gravity force is 
neglected and the vehicle is allowed to heave and rotate freely. The second case is that of a free flight simulation. 
Mechanisms are proposed to explain the stability behavior of the vehicle. 

II. Modeling 


A. Vehicle Definition 

Figure 1 shows an aeroshell with an isotensoid attached IAD during atmospheric descent. Because the altbody 
of the aeroshell is mostly enclosed by the IAD, only the geometry of the aeroshell forebody, a 70-degree sphere 
cone, is considered in this analysis for the purposes of aerodynamics. The IAD is a scaled version of an isotensoid 
shape analyzed and tested by Mikulas and Bo lion. ; The shape is axisymmetric and does not include ram air inlets or 
a burble fence. The IAD mass properties were determined using sizing relationships presented by Clark et al. 4 and 
the mass moment of inertia properties were estimated. Table 1 contains geometric, mass, and mass moment of 
inertia properties for both bodies with the geometry definition given in Figure 2. 

In order to add a rotational degree of freedom between the aeroshell and the IAD, they are treated as separate 
bodies that are attached to each other at a single point. Both bodies remain independently rigid and therefore 
compatibility between the perimeter of the aeroshell and the IAD is not enforced. The interface is modeled as a 
spring-damper which constrains the relative motion of the two bodies. Figure 3 illustrates the flexible configuration 
at an arbitrary pitch difference between the two bodies. 

B. Equations of Motion 

The equations of motion of the vehicle are derived in the body-fixed coordinate frames attached to each body 
with the origins coincident at the attachment point. For each body, the pitch angles are defined as positive above the 
local horizon (see Fig. 1). Figures 4 and 5 illustrate the forces and moments acting on each individual body at the 
attachment point as well as the definition of both body axes of the system. 

Assumptions made in deriving the equations of motion are as follows: 

1) The system has four independent degrees of freedom: u, w, 9 , and A 0 . 

2) The x-z and x ’-z ’ planes are coplanar and all motion is constrained to these planes. 

3) The center of mass of the aeroshell is fixed but may be offset from the aeroshell’s centerline. 

4) The center of mass of the IAD is fixed on the IAD’s centerline. 

5) The attachment between the IAD and the aeroshell is modeled as a single point (i.e., points C and C’ in Figs. 

4 and 5 are coincident). 

6) The planet is flat. 

7) There is no wind. 

8) The acceleration of gravity, g, is constant. 

The translational equations of motion as derived from Newton’s second law at a point other than the center of 
mass are 5 


■ F a + R x - m S sin<9 = m{u + Qw + Qz CM - Q 2 x CM ) 

(1) 

F n + R z + mg cos 0 = m(w -Qii- Qx CM - Q 2 z CM ) 

(2) 

- F a + R y - m gs\nO = m (w +Qw -Q 2 x CM ) 

(3) 

-F n + R, + m gcos0 = m (w -Qu -Q x CM ) 

(4) 


It is desirable to eliminate the interface forces from the translational equations of motion because they are forces 
internal to the aeroshell/lAD system. This can be accomplished by recognizing that 
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R x = —R x cos(<9 — 0 )-R : sin (# -6 ) 


(5) 


R- = +R-, sinl0 - 


in [o -0 )- R z cos(# - 0 ) 


(6) 


Substituting Eqs. (5) and (6) into Eqs. (3) and (4) allows the interface force to be removed from the translational 
equations of motion by solving the resulting system of equations. Rearranging the resulting equations yields the 
translational equations of motion for the aeroshell as 


\-F' A -[F A +mgsmO + m(gw + Qz CM -Q 2 x CM \ cos^-tf')] 


m cos 




■ [f v - mg cos 0 + m (w - Qu - Qx ( 
-mgsmO -m ( u +Qw -Q 2 x CM J 


K CM 


Q Z CM 


)]sin(6»-6») 


m cos 6* - 




-F, 


+ F 4 + mg sin 0 + m(u + Qw + Q z cm 


- J'jr T I' A X SHIP -r m\u t v t^ cm -Q Xq U jjsiiqo 

- [f w - mg cos 0 + m(- Qu - Qx CM -Q 2 z CM )]cos(^ - 0 ) 


j sin(^ - 6 )] 


+ m g cos 0 - m (w —Q u — 


‘ Q x cm ) 


(7) 


( 8 ) 


The two equations for it and w can be obtained by recognizing that at the attachment point between the aeroshell 
and the IAD 


v c = v c . 


— =—\ r 


dt 


dt 


(9) 

( 10 ) 


By carrying out the derivatives on both sides of Eq. (10) and expressing the aeroshell body axis acceleration in the 
IAD coordinate frame, the translational equations of motion for the IAD are 


u = (u + 0w)cos(t? -6 )+ (w — Qu) sin(# -6 )-Qw 

w = —(u + 0vy)sin [d-0 )+ (w - Qu)cos[p - 0 )+Q u 
where the velocity of the IAD, written in terms of the velocity of the aeroshell, is 

u = u cos(# - 9 )+ u’sin(# - 6 ) 
w = -u sin(# -0 )+ wcos(# - 6 ) 


( 11 ) 

( 12 ) 

(13) 

(14) 


Equations (11) and (12) can be substituted into Eqs. (7) and (8), yielding the two degree of freedom equations of 
motion in translation in the aeroshell body-fixed frame. 

The rotational equations of motion for the aeroshell and the IAD are taken about the coincident origins of the 
body-fixed frames and are given by 


m r + M aero,c ~ m g{ x CM cos 6» + z CM sin 0) = o[l ^ + m{z 


CM “ r A CM 


2 + Xru )l - + mz CMU + m Q{ wz CM +ux Cm) ( 15 ) 


-CM U 


~M R + M aeroC - mgx CM cos 6 =Q [l\ y + m x'c M \~ m x CM w + m Q ux' CM 


(16) 


By rearranging Eqs. (15) and (16), the rotational rate equations of motion for the aeroshell and the IAD are given by 
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Mr is the interface moment between the aeroshell and 1AD. For this derivation the interface moment will be 
modeled as a torsion spring and damper of the form 

M R =-K(e-e')-?{Q-Q) (19) 

The explicit form of the equations of motion given by Eqs. (7), (8), (17), and (18) are obtained by casting the 
equations in matrix form and solving for the state variable time derivatives. Doing so yields the equations of motion 
as a function of the state variables u, w, Q , Q , 8 , and 0 . The equations of motion for the pitch angle of both 
bodies are 


0 = Q 

(20) 

0=Q 

(21) 

Because the calculation of atmospheric density requires the altitude of the vehicle as input, another equation of 
motion is required to obtain the vertical location of the vehicle from the surface of the planet in inertial space. This 
equation can be obtained by resolving the body-fixed velocity components of the aeroshell vertically using the 
equation 

h = usind- wcos8 

(22) 

C. Atmospheric Model 

The calculation of density is accomplished using an exponential model which is a curve 
descent of the Viking Lander. 6 The model has the general form 

fit of data from the 

p(h)= A>expf--^rj 

(23) 


where in the Martian case the reference density, p Q , is 0.0216 kg/m 3 and the scale height, H, is 9.76 km. 

D. Aerodynamics 

The aerodynamics of both the aeroshell and the IAD were calculated using modified Newtonian aerodynamics. 
While this method of analysis is typically reserved for the hypersonic regime of flight, the aerodynamics of an 
isotensoid IAD remain roughly constant through the supersonic regime as it approaches the hypersonic limit. 1 
Therefore it is assumed that this method for determining the aerodynamics of the IAD remains valid in the 
supersonic regime spanning Mach 2.5 to 5. Comparison of modified Newtonian aerodynamic results for the 
isotensoid shape with experimental data 3 verifies the validity of this assumption. 

Aerodynamic moments for both bodies are calculated about the attachment point C/C’ . The reference area and 
length used are the projected area and diameter, respectively, for each body (see Figure 2). Note that the calculated 
aerodynamic forces and moments for the IAD do not include the depicted flat front panel, but the only external 
curved boundary. The aerodynamic moment coefficient trend with respect to angle of attack for both bodies 
indicates they are both statically stable about the point C/C’ (see Table 2). The coefficient of both moment and 
normal forces with respect to angle of attack are linear within the range of interest. Axial forces are a function of 
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angle of attack for both bodies with the maximum value at zero angle of attack indicated in Table 2. The value for 
base pressure coefficient was set to zero during the aerodynamic analysis. Equation (24) expresses the calculation 
of total aerodynamic moment due to static and dynamic aerodynamics. 


M 


aero 



a + C m 


Qd_ 
2 2V 



(24) 


The effect of pitch damping derivatives (i.e., C m ^ ) are included in this study. For this analysis, pitch damping is 

assumed to be constant (i.e., not varying with angle of attack) for each body. The stable (i.e., negative) value of 
pitch damping is obtained using modified Newtonian aerodynamics. For the unstable value of pitch damping, the 
sign of the stable result was changed from negative to positive. Neutral (i.e., zero) pitch damping is used by default 
in this analysis with the effects of stable and unstable pitch damping determined in a parametric study for the free 
flight case. 

E. Non-dimensionalization of Interface Quantities 

Non-dimensionalization of the interface stiffness and damping makes it easier to interpret the results. The 
interface stiffness is non-dimensionalized using static aerodynamics while the interface damping is 
non-dimensionalized according to the damping ratio found in literature. 6 7 The non-dimensional interface stiffness is 
given by the equation 


C ma q xfi (S' +S)d' 

where the denominator represents an effective stiffness due to aerodynamic forces. The aerodynamic moment 
coefficient is that of the combined bodies about the attachment point. The dynamic pressure is the initial value for 
the problem being investigated. The reference area and length are the projected area and diameter of the entire 
vehicle. Therefore, Eq. (25) provides a transformation from the dimensional spring stiffness to a 
non-dimensionalized stiffness for a given geometry and initial flight condition. 

Non-dimensionalizing the interface damping is achieved by using the damping ratio, defined as 


v = 


n 

2 yj ld eff 


(26) 


where the denominator is the critical damping for a torsion spring. 8 The effective moment of inertia is related to the 
relative oscillation mode between the aeroshell and the IAD as they oscillation in a vacuum. The effective moment 
of inertia is defined to be 


: eff ~ 


1 -- 


yy 


I yy +I yy J 


(27) 


Substituting Eq. (27) into Eq. (26) yields the non-dimensional spring damping. Note that the presence of the spring 
stiffness in the non-dimensional damping term causes a correlation between the selections of the two 
non-dimensional quantities. This correlation could be avoided by non-dimensionalizing the damping against a 
relevant non -zero value of C m ^ (e.g. its value at zero angle of attack). 


III. Results 


A. Simulation Validation 

In order to compare the behavior of the multibody free flight simulation to already-existing solutions, the 
non-dimensional interface stiffness and damping were set arbitrarily high to cause the vehicle to behave as a single 
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rigid body. The numerically integrated angle of attack history of the vehicle was compared to a closed form 
solution’ for a Mars flight case. Both analyses yielded identical results. 

B. Constant Velocity Results 

In this modification to the free flight case, a virtual thrust is applied at the attachment point to counteract the drag 
experienced by both bodies at a particular instant. By adding this virtual thrust, the vehicle is allowed to have a 
steady airspeed after the initial transients are passed. Because a virtual thrust is used instead of a virtual sting, the 
vehicle is still allowed to heave. The effect of gravity is also removed from the simulation leaving only the effects 
of the interface and aerodynamics. The initial conditions for the simulation are given in Table 3 with the density 
being held constant. Each case was run for five seconds. 

Stability regimes for the constant velocity simulation are obtained by isolating the longest period motion of the 
angle of attack of the aeroshell with respect to time. For each run representing a unique combination of k and rj , a 
curve fit is applied to the numerical results of the aeroshell angle of attack time history. Because for some interface 
designs there is a secondary, short-period oscillation, the nonlinear curve fit must have the form 

a(t) = A x e Blt cos^t + <f\ )+ A 2 e B - cos (co 2 t + (/> 2 ) (28) 

The exponential decay constant corresponding to the longer period motion (i.e., the w l or co 2 with lower value) is 
used in the constant velocity case as a measure of stability. Acceptable fits were determined to be those with a 
coefficient of determination to the numerical results in excess of 0.99. 

Figure 6 shows contours of the decay constant over a large range of non-dimensional interface stiffness and 
damping values for C m ^ =0. Two major regions of stability associated with the constant velocity case can be 

discerned. While both regions are characterized by statically stable behavior, the region spanning high 
non-dimensional interface stiffness values exhibits dynamic instability of the aeroshell angle of attack as time 
progresses (i.e., B > 0). In the other region, the angle of attack amplitude of the aeroshell is dynamically stable (i.e., 
B < 0) and therefore decays as time progresses. As designs move toward the interior of this region, B becomes 
increasingly negative and thus the time history of aeroshell angle of attack becomes more highly damped. 

C. Free Flight Results 

The free flight simulation reflects the behavior of the vehicle as it decelerates in the Martian atmosphere while 
under the influence of the planet’s gravity field. The initial conditions for this case are given in Table 4 and the 
simulation was terminated once Mach 2.5 was reached or the angle of attack of either the aeroshell or the 1AD 
exceeded 30 degrees. 

Figure 7a qualitatively shows three distinct regions of stability associated with the free flight case for C m<) = 0. 

The region that spans high values of non-dimensional interface stiffness is characterized by statically stable, 
dynamically unstable oscillation of the aeroshell angle of attack. This is the same behavior as observed in the 
constant velocity simulation at high non-dimensional interface stiffness values. In the region spanning moderate 
non-dimensional interface stiffness values there is statically stable oscillation of aeroshell angle of attack which 
decays to a lower amplitude with respect to the beginning of the simulation. This behavior is defined as being 
statically and dynamically stable. In the region encompassing low non-dimensional interface stiffness values, the 
aeroshell is initially driven to angles of attack exceeding 30 degrees. For this analysis such behavior is defined as 
being divergent, even though for some cases the aeroshell angle of attack may return to lower values due to the 
restoring influence of aerodynamic and interface moments when the system oscillations grow to large amplitudes. 

The effect of gravity is negligible in the behavior of the system as removing gravity did not affect the stability 
boundaries seen in Fig. 7a. Applying a stable pitch damping ( C m<) < 0 ) to both bodies had the effect of removing 

the boundary separating the dynamically stable and unstable regions for the ranges of k and fj considered and 
replacing it with a single region of static and dynamic stability as seen in Fig. 7b. As non-dimensional interface 
stiffness is decreased the aeroshell stability is dominated by the effect of the drag of the 1AD. As such, applying a 
stable pitch damping did not affect the boundary between the statically and dynamically stable region and the 
divergent region. Applying an unstable pitch damping ( C m > 0 ) to both bodies had the effect of lowering the 

boundary between the statically stable but dynamically unstable and the statically and dynamically stable regions as 
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seen in Fig. 7c. Again, the boundary between the statically and dynamically stable region and the divergent region 
remains unaffected due to the dominance of the IAD drag at low values of non-dimensional interface stiffness. 


D. Stability Mechanisms 

All cases discussed in this section have neutral pitch damping ( C m =0). With high values of the non- 

dimensional interface stiffness the aeroshell and IAD exhibit statically stable but dynamically unstable trends of 
aeroshell and IAD angle of attack as shown in Fig. 8. This behavior is observed both in the constant velocity and 
free flight simulations. For the constant velocity simulation, the generation of lift causes dynamically unstable 
behavior due to the heaving the vehicle. For the free flight case, both the generation of lift and the decreasing 
dynamic pressure cause dynamically unstable behavior. 2 

Reducing the non-dimensional interface stiffness has the effect of creating statically and dynamically stable 
angle of attack time histories for the aeroshell and IAD. In this region, the vehicle interface is a dominant factor in 
the stability of the vehicle. Figure 9 shows such a case from the free flight simulation. The difference in 
aerodynamic moment between the aeroshell and the IAD induces an initial difference in pitch and pitch rate between 
the bodies. The lag in pitch rate activates the damping mode of the interface causing dissipation of rotational energy 
from the system. This dissipative effect can be removed if non-dimensional interface stiffness is held constant and 
the non-dimensional interface damping is reduced significantly, as shown in Fig. 10. Without the effect of interface 
damping, the same dynamically destabilizing effects as observed in the highly stiff limit are once again dominant. 

The region of divergent behavior is observed only in the free flight case where both the non-dimensional 
interface stiffness and damping have low values. Figure 11 shows a case in this region from the free flight 
simulation. When the non-dimensional interface stiffness and damping are in this region, the two bodies are given 
significant freedom to initially rotate independently of each other. Because the aerodynamic moment acting on the 
IAD is much greater than that acting on the aeroshell, the IAD quickly rotates into the oncoming flow and oscillates 
about an angle of attack of zero while the aeroshell remains at a positive angle of attack. The drag force of the IAD 
is transmitted to the aeroshell through the interface with an angular bias with respect to the aeroshell centerline. As 
the aeroshell center of mass is located on the centerline behind the attachment point, a pitch-up moment is 
experienced by the aeroshell. Figure 12a shows the components of the moments on the aeroshell about its center of 
mass for the case in Fig. 11. Because the moment due to the drag of the IAD becomes stronger than the aeroshell’s 
combined aerodynamic and reaction moments, the aeroshell angle of attack diverges. This can be verified by Fig. 
12b, which gives the sum of all the moments in Fig 12a. Note that this behavior is not witnessed in the constant 
velocity simulation due to the addition of the drag-countering virtual thrusts to each body. 

IV. Conclusions 

This paper has presented the equations of motion for a four degree of freedom vehicle consisting of an aeroshell 
with an attached IAD. Using these equations, along with aerodynamic data generated by numerical computation, the 
effect of varying the non-dimensional interface stiffness and damping on vehicle stability was determined. It was 
shown that, depending on the combination of non-dimensional interface stiffness and damping, the effects on the 
static and dynamic stability of the aeroshell can be significant. For high values of the non-dimensional interface 
stiffness the system behaves as a single rigid body and exhibits the characteristic statically stable but dynamically 
unstable time history of the angle of attack of rigid bodies. For moderate values of non-dimensional interface 
stiffness, the system exhibits statically and dynamically stable behavior of aeroshell angle of attack due to the 
dissipating effect of the interface damping. For low values of interface stiffness, the system exhibits divergent 
behavior of aeroshell angle of attack due to the biased pointing of the IAD drag force. 

While no provision has been made in this study to determine physically realistic values of non-dimensional 
interface stiffness and damping, the framework developed is able to accommodate such data. If provided with 
specific values of stiffness and damping, the simulation will predict the stability behavior of the vehicle during 
flight. Therefore, when combined with a fluid-structure interaction code, this simulation framework is a useful tool 
for evaluating the stability of IAD designs. 
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Table 1. Mass and geometry properties of the aeroshell and 1AD. Coordinate systems for center of mass location are 

defined in Figures 4 and 5. 


Parameter 

Aeroshell 

IAD 

m 

2648 kg 

29 kg 

lyy 

2863 kg-m 2 

294 kg-m 2 

d 

4.5 m 

10.18 m 

s 

15.9 m 2 

65.5 m 2 

*CM 

-0.44 m 

-2.66 m 

ZCM 

0.00 m 

Not Applicable 


Table 2. P redicted aerodynamic quantities for the aeroshel l and IAD. 


Parameter 

Aeroshell 

IAD 

c 

-0.2007 rad 1 

-0.3881 rad 1 

c N 

a 

0.2761 rad' 1 

0.7786 rad 1 


f -0.360 rad 1 

-0.311 rad 1 

C„. H 

0 

0 

rn Q 

[_ 0.360 rad' 1 

0.311 rad 1 

Max Cj 

1.5498 

0.7675 


Table 3. Initial conditions for constant velocity simulation. 

Parameter 

Value 

M 

4.5 

a, a 

5 deg 

e 

0 deg 

ao 

0 deg 

p 

2.78><10" 3 kg/m 3 

Q.Q’ 

0 deg/sec 

Table 4. Initial conditions for free flight simulation. 

Parameter 

Value 

M 

5 

a, a 

5 deg 

9 

0 deg 

a e 

0 deg 

h 

20 km 

Q, Q’ 

0 deg/sec 
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Planet Surface 


Figure 1 . Aeroshell with isotensoid 1AD in flight. 


Profile View 


Projected View 


Figure 2. Geometry definition for the aeroshell and IAD 


American Institute of Aeronautics and Astronautics 




Figure 3. Flexible aeroshell and 1AD system. Moments due to the spring 
and damper at the interface point depend on the relative angle and its time 
derivative, respectively. 



Figure 4. Aeroshell with external forces, interface forces and 
moments, and coordinate system definition. 
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Figure 5. IAD with external forces, interface forces and moments, 
and coordinate system definition. 



Statically Stable 
Dynamically Unstable 
(B > 0) 


Statically Stable 
Dynamically Stable 
(B < 0) 


n 


Figure 6. Contours of the decay constant, B, are shown as a function of non-dimensional 
interface stiffness and damping with C m ^ = 0 . 
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71 


Figure 7a. Qualitative stability regimes for aeroshell angle of attack in free 
flight with C m<) = 0 . 



71 


Figure 7b. Qualitative stability regimes for aeroshell angle of attack in 
free flight with C mg < 0 . 
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1 


Figure 7c. Qualitative stability regimes for aeroshell angle of attack in 
free flight with C„ !y > 0 . 



Figure 8. Statically stable but dynamically unstable angle of attack 
oscillation trends for free flight case with k = 1000, fj = 2, C m ^ =0. 

Note that the time history of angle of attack for the aeroshell and 1AD are 
the same and that Ad remains at zero. 
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Time (sec) 

Figure 9. Statically and dynamically stable angle of attack oscillation 
trends for free flight case with k = 3, r) = 2, C„, o = 0 . 



Figure 10. Statically stable but dynamically unstable angle of attack 
oscillation trends for free flight case with k =3, r) = 0.05, C m ^ = 0 . 
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Figure 11. Divergent aeroshell angle of attack oscillation trend for free 
flight case with k =0.001, rj = 2, C m =0. 



Figure 12a. Components of the moment about the aeroshell center of 
mass for a free flight case with k = 0.001, fj = 2, C =0. The 

moment due to gravity is neglected. 
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Figure 12b. Total moment about the aeroshell center of mass for a free flight 
case with k = 0.001, rj = 2, C m ^ =0. The moment due to gravity is 

neglected. 
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